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Abstract. We consider a physical scenario for ultra-relativistic heavy-ion 
collisions where, at the early stage, only transverse degrees of freedom of partons 
are thermalized, while the longitudinal motion is described by free streaming. 
When the energy density of the partonic system drops to a certain critical value, 
the partons hadronize and the newly formed hadronic system freezes out. This 
sudden change is described with the help of the Landau matching conditions 
followed by the simulations done with THERMINATOR. The proposed scenario 
reproduces well the transverse-momentum spectra, the elliptic flow coefficient 
V2, and the HBT radii of pions and kaons studied at RHIC (Au+Au collisions 
at y/SNN = 200 GeV). It also reproduces quite well the transverse-momentum 
spectra of hyperons. 



PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 

The spacetime evolution of matter produced in the uhra-relativistic heavy-ion 
coUisions is described most successfuUy by the perfect-fluid hydrodynamics [H 121 13] ■ 
With the initial conditions inferred from the Glauber model, the hydrodynamic 
approach describes very well the hadron transverse-momentum spectra and the elliptic 
flow coefficient V2. With the Gaussian initial conditions, the hydrodynamic approach 
describes also the HBT radii [4]. 

Unfortunately, the use of the hydrodynamics is still challenged by the problem of 
early thermalization ~ the correct description of the data requires an early initialization 
time for hydrodynamics 5j, which implicitly assumes a very early equilibration time. A 
possible solution to this problem was suggested in Ref. JHJ , where the idea was proposed 
that only transverse degrees of freedom are thermalized, while the longitudinal motion 
is described by free streaming (a physical picture explained naturally in the string 
models [TJ [5]). In this framework, called below the transverse hydrodynamics, one 
can obtain the parton transverse-momentum spectra and V2 which are consistent 
with the data [SI [9] [§|. In this paper we generalize this approach by inclusion of 
the sudden-hadronization transition modeled with the help of the Monte Carlo code 
THERMINATOR 12J. 

X Supported in part by the Polish Ministry of Science and Higher Education, grant N202 034 32/0918. 
§ We note that our formulation differs from the formalism of the transverse hydrodynamics introduced 
originally in [101 111) . For a more detailed discussion of this issue see [6]. 
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The equations of the transverse hydrodynamics foUow from the energy-momentum 
conservation law, (9^r|"' = 0, with the energy-momentum tensor [13] 

= — [(£2 + P2) U^U-" - P2 (g'"- + V^V) ] , (1) 

T 

where r = \/W^^z^ is the longitudinal proper time and no describes the density of 
transverse clusters in rapidity. The clusters are identified with the groups of partons 
having the same rapidity. They are two-dimensional (2D) objects described by the 
2D thermodynamic variables: £2, S2 and T2 (2D energy density, pressure, entropy 
density, and temperature, respectively). The 2D thermodynamic variables satisfy the 
standard thermodynamic identities: £2+^2 = T2S2, ds2 = T2ds2, and dP2 = S2dT2 
(the baryon chemical potential is set equal to zero). The energy- momentum tensor 
([1]) includes also the two four- vectors, 

C/^ = (mo cosh?7, Uy, "0 sinh?]), = (sinh?7, 0, 0, cosh?]), (2) 

where = {vP^u^^Q) is the hydrodynamic flow in the plane z = 0, while 
ri = 1 / 2hi{{t + z) / {t — z)) is the spacetime rapidity. The four-vectors C/^ and 
satisfy the following normalization conditions: U^U^ = 1, V^V^ = — 1, 
U^Vfj^ = 0. The four-vector corresponds to the flow four-velocity in the standard 
hydrodynamics. On the other hand, the term V^V^ in ^ is responsible for vanishing 
of the longitudinal pressure, i.e., in the local rest-frame of the fluid element, where 
U^' = (1, 0, 0, 0) and V^' ^ (0, 0, 0, 1), we have 



/ £2 \ 

no F2 

T F2 

\ y 



(3) 



It should be emphasized that such a structure of the energy-momentum tensor 
appears in the theory of the color glass condensate and glasma for t ^ ^/Qs, where 
Qs ^ I GeV is the saturation scale expected at RHIC [TU [15] . This gives support for 
consideration of the transverse hydrodynamics as the appropriate description of the 
early evolution of matter. 

The equations of the transverse hydrodynamics are solved numerically with the 
equation of state P2 = £2/2 = i^gT^ C(3)/(27r), where i^g = 16 [13]. Our analysis is 
restricted to the central rapidity region (z w r/ w 0) , where the created system may be 
treated as boost- invariant. The initial conditions assume that the 2D initial energy 
density at t = Tj = 1 fm is proportional to the mixture of the wounded-nucleon 
density, pw ix±), and the binary-collision density, ps (a;±), 

£2 (Ti, x±) = 9 2V. ysw ^ —^p^ (f^) -I- npB {x±) , (4) 

TT Z 

where, following the PHOBOS studies of the centrality dependence of the hadron 
production [16j, we take k = 0.14. The normalization constant required in (|4]) 
determines the 2D initial central temperature of the system, T2i ~ T2{Ti, 0). 
The Landau matching conditions have the form 

TrC/. = T^'^U,, (5) 
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where Tg is the standard energy-momentum tensor of the perfect-fluid hydrodynam- 
ics 

T^-" = (£3 + PsWU'' - P^gf^r (6) 

Here £3 and P3 are the three-dimensional (3D) energy density and pressure. Equations 
© and © give 

no 



T 



-£2 = £3, (7) 



which should be supplemented by the requirement of the entropy growth, 

— S2 < S3, (8) 

r 

where S3 is the 3D entropy density. Dividing both sides of Eqs. ^ and ([8]) one obtains 

> (9) 
2 S3 

In our calculation, £3, P3, and S3 are interpreted as the quantities characterizing 
the hadron resonance gas in equilibrium at the 3D temperature T^f. This means 
that the use of the Landau matching conditions is a formal scheme to describe the 
fast hadronization process; if the energy density of the system drops to £3(T3f), the 
transverse gluonic clusters change into a locally isotropic resonance gas. 

The considered model has three free parameters: T2i, T^f, and uq (the initial 
time T[ is always set equal to 1 fm). In the fitting procedure we first assume uq — 1 
and try to adjust and T^i in such a way that a good description of the data 
is achieved. In this step, the energy-conservation condition ([7]) is used as the only 
constraint. In the next step we rescale the 2D temperature in such a way that the 
equality T2 = 3 £3/(2 S3) holds at the initial time r = Tj. In this way we make certain 
that the condition ([9]) is fulfilled on the whole transition hypersurface. The rescaling 
of the temperature changes the energy density on the 2D side, however, this may be 
compensated by the appropriate change of no (the flow profile is not affected by the 
rescaling of the temperature). In this way we are able to satisfy finally the conditions 
© and dHl). 

Our description of the hadronization process is quite radical. Yet, we add to 
it another radical assumption: the evolution of the isotropic (locally equilibrated) 
hadronic system is so short that the hadronization process may be followed 
immediately by freeze-out. Technically, this assumption is realized by the Monte-Carlo 
simulations done with THERMINATOR, where the 3D hadronization temperature 
T3f is identified with the freeze-out temperature and the hadronization hypersurface 
becomes equivalent to the freeze-out hypersurface. 

Certainly, our method used to switch from the purely transverse hydrodynamic 
expansion to the locally isotropic system of hadrons and further to freeze-out is 
very much simplified. In particular, one expects that the isotropization process 
has a gradual character. More elaborate approaches would describe this kind of 
transformation using kinetic theory or dissipative hydrodynamics, for example, see 
[17l [iHl [19] . Alternatively, one may think of several direct developments of our scheme, 
for example, the transition criterion may be changed and/or the time evolution of the 
isotropic phase may be extended. Such developments are the subject of our current 
studies. Anyway, in our opinion such more realistic smooth changes may be effectively 
replaced by the sudden but delayed processes. 





Figure 2. Elliptic flow coefficient of pions and kaons (solid line) compared 
to the experimental PHENIX data (dashes) 



The results of our numerical calculations with the initial profiles corresponding to 
the non-central (c — 20-30%) Au-|-Au coUisions at the beam energy ^snn = 200 GcV 
are shown in Figs. [T]-|4l The values of the parameters used in the calculations are: 
T2i = 590 MeV, Tgf = 165 MeV, and no = 0.43 [jj] . We observe that our model 
reproduces well the transverse- momentum spectra of pions, kaons, and protons, Fig.[l] 
and also the pion-|-kaon elliptic flow, Fig. [21 both measured in Au-t-Au collisions by 
PHENIX [5T1 [25] . The values of A indicate the average relative difference between 
the experimental data and the model results. 

In Fig. [3] we show the model pionic HBT radii (solid lines) compared with the 
experimental STAR data (dots) [53]. The radii were calculated with the help of 

II These parameters are slightly different from those used in 1201 . where we did not study the 
production of strangeness and the equation of state was of the form P2 = £2/2 = VgT^/{2n). 
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Figure 3. Model results for the pionic HBT radii (solid lines) compared 
to the experimental STAR data [23] (dots). The radii are plotted as the 
functions of the average transverse momentum of the pion pair kr- The 
model parameters are the same as in Figs. [T]and[2] 




Figure 4. Model results for hyperons compared to the experimental STAR 
data [231 (dots). The model parameters are the same as in Figs. [T]and[2] 
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THERMINATOR which uses a two particle method with Coulomb corrections |25j . 
The overall agreement with the data is quite good. The theoretical radii are by 
about 20% smaller than the experimental values, but their dependence is very well 
reproduced. Moreover, the model ratio i?out/-Rsidc matches very well the experimental 
result. Slightly smaller theoretical values indicate that the real isotropization and 
freeze-out processes take longer time. This is, of course, an expected feature. 

The high freeze-out temperature obtained in our model suggests that one may 
be able to reproduce correctly the spectra of hyperons, as in the case of the single- 
freeze-out model [55]. The hyperons, due to their small interaction cross sections, 
are believed to decouple rather early at the temperature close to the typical chemical 
freeze-out temperature. Indeed, our calculations presented in Fig. [4] show that the 
hyperon transverse-momentum spectra are reproduced reasonably well. 

In the conclusions we state that the transverse-hydrodynamics model 
supplemented with the sudden isotropization transition and freeze-out is able to 
describe uniformly the soft pion and kaon data (transverse-momentum spectra, U2, 
and the HBT radii). This brings further evidence that the assumption of the very 
fast 3D equilibration of matter produced in the relativistic heavy-ion collisions may 
be relaxed — an earlier study showed that the RHIC data may be explained by 
the model where the standard hydrodynamic phase is preceded by the parton free 
streaming [57]. Our analysis of the HBT radii indicates, however, that the ultimate 
3D equilibration is necessary to reproduce the experimentally measured value of i?iong- 
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